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Abstract. We report new results on the cluster of galaxies Abell 970 obtained from X-ray observation with 
BeppoSAX. Our analysis of the BeppoSAX MEGS and LEGS data in the range [0.15-10] keV reveals a mean 
cluster gas temperature of kT = 4.461q']^ 5, a metallicity of Z = O.SIIq q^Z©, and an interstellar hydrogen 
absorption density column of A^h = O.OSlg gY x 10^°cm“^. Moreover, we obtained azimuthally averaged radial 
profiles of these quantities. Our results are consistent with the hypothesis that Abell 970 has been disturbed by 
a past merger or by the ongoing merger process of a substructure, that put the cluster out of equilibrium. This 
is also demonstrated by the offset between the gas and galaxy distributions. Gombining the X-ray data with a 
recently published analysis of new galaxy radial velocities, we conclude that a subcluster 8 aremin to the NW is 
falling into Abell 970 and will merge in a few Gigayears, thus disturbing Abell 970’s newly acquired equilibrium. 
The high a-elements/iron ratio that we derive for this cluster supports the hypothesis of early intracluster medium 
enrichment by Type II supernovae. 
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1. Introduction 


Rich clusters of galaxies are the last structures to col¬ 
lapse in the Universe. When isolated and within its virial 
radius, the cluster contents should already be in quasi¬ 
stationary equilibrium by now. However, according to the 
hierarchical structure formation scenario, rich clusters are 
formed by the accretion of dark matter dominated haloes 


(White & Reeg 1 

978 

and, e.g., the W-body simulations of 

[Evrard et al. 

200 

2). These recurring merger events will af- 


feet the cluster, moving it away from an equilibrium state. 
Therefore, the study of the morphology and dynamics of 
cluster of galaxies can give valuable information on struc¬ 
ture formation in the Universe. 


The hot {kT « 2-10 keV) and tenuous (central density, 
riQ 10“^cm“^) X-ray emitting gas found in rich clusters 
of galaxies is an excellent tool to probe the cluster proper¬ 
ties and past histor y. As it was alre ady recognized almost 
three decades ago ( Lea et al. 1973 ), the short relaxation 
time of the gas, and the fact that the main emission mech¬ 
anism is the optically thin thermal bremsstrahlung, make 
it an excellent tracer of the cluster gravitational potential 
well. 

Besides the temperature and density, the X-ray ob¬ 
servation of the intracluster medium (ICM) can be used 
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to measure the metal abundance (hereafter metallicity) of 
mainly iron, but also of some a-elements. Such observa¬ 
tions can contribute to our understanding of gas enrich¬ 
ment in clusters, which depends on the early evolution of 
galaxies in clusters and their interaction with the intra¬ 
cluster medium. 


Here, we present an X-ray study of Abell 970 {z = 
0.0587, richness clas s J? = 0, B-M type HI) following the 
work of ^odre et al. ( 2001 , Paper I hereafter). 

This cluster is in a rather dense environment, being 
part of the Sextans supercluster (number 88 in the cata¬ 


logue of Einasto et al. 1997, and number 378 in the cat¬ 


alogue of Kalinkov et al. I998| , which also includes Abell 
978, 979, and 993). It has moderate X-ray luminosity and 
it was observ ed with the Einstein I maging Proportional 
Counter (e.g. jjones fc Forman 1999 ); it w as also detected 
in th e ROSAT Bright Survey Catalogue ( ^chwope et al ' 


^2000 ), although a pointed observation has never be e n car- 
ried out with ROSAT. In the study of White et al. ( 1997 ) 
using Einstein IPC data, it has been found that Abell 970 
has at most a weak coo ling flow, of about 20Mq yr“^ (see 
also Loken et al. 1999). 


In the dynamical analysis made in Paper I, we found 
that the galaxies in Abell 970 have a mean velocity dis¬ 
persion of 845 km s“^ (and up to ^ 1000 km s“^ in the 
centre). This analysis also indicated that this cluster must 
have substructure and is out of dynamical equilibrium. 
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This conclusion is also supported by an offset in the peaks 
of the surface density distribution and the X-ray emission 
obtained from an Einstein IPC image. Further support 
came from a large-scale velocity gradient in the cluster. We 
also found a discrepancy between the dynamical masses in¬ 
ferred from the virial theorem and those inferred with the 
X-ray emission, which is expected when the galaxies and 
the gas inside the cluster are not in virial (or hydrostatic) 
equilibrium. Here we present new results on the distribu¬ 
tion of temperature, density, and abundances of the X-ray 
emitting gas, as well as its a-elements/iron ratio. 

This paper is organized as follows. We describe in 
Sect. H the data obtained with BeppoSAX. In Sect. || we 
present the X-ray fluxes and luminosities and determine 
and analyse the temperature and metallicity radial pro¬ 
files. The cluster gas and dynamical masses are estimated 
in Sect. In Sect. || we discuss the possible merging sce¬ 
nario of Abell 970, and in Sect. || we give the individual 
metal abundances and discuss the consequences for the 
ICM enrichment modelling. 


2. The Data 


Abell 970 was observed in July 2000 by the BeppoSAX 



33 619 s for the MEGS and LEGS, respectively. 

We restricted the scientific analysis to the 8 arcmin 
(~ 760 h^Q kpc at the distance of Abell 970) circle cen¬ 
tered at the aim point of the LEGS and MEGS instruments 


in order to prevent severe vignetting effects (Gusumano & 
Mineo| 1998). Moreover, the MEGS entrance window is 
supported by a beryllium “strongback” that is opaque to 
low energy photons {E 5 keV) and is not taken into ac¬ 
count in the effective area files (Fiore, private communica¬ 
tion). The strongback structure affects detector positions 
between 7.9' to 13.3' ( pusumano fc Minec 1998). Since 
we are analysing data inside the MEGS window struc¬ 
ture, we opted for combining both MEGS 2 and 3 event 
files in order to improve the signal-to-noise ratio, instead 


of ana lysing them separately as De Grand! & Molendi 
( |2002 ), which would be essential for an analysis extending 
beyond ~ 8'. 

Figure shows the X-ray isocontours superposed on 
a “Second Epoch Survey” DSS optical imageQ. The event 
file coordinates were corrected for the systematic error re¬ 
ported by the BeppoSAX team.^. After correction, we ver¬ 
ified a very good matching of both the MEGS and LEGS 


^ The “Second Epoch Survey” of the southern sky was made 
by the Anglo-Australian Observatory (AAO) with the UK 
Schmidt Telescope. Plates from this survey have been digitized 
and compressed by the STScI UK Schmidt Telescope and the 
Digitized Sky Survey. 

^ cf. www.asdc.asi.it/bebbosax/coord_correction.html. 


images to the Einstein IPG image (see Paper I), includ¬ 
ing the small X-ray peak at the SE, which, as pointed 
out in Paper I, is associated to a pair of interacting clus¬ 
ter galaxies (numbers 37 and 38 of Paper I catalogue: 
a = 10hl7“37", S = -10°46'02", J2000). The distance 
between the optical position of the pair and the associ¬ 
ated X-ray peak is less than 15". 

A comparison of Figure with the Einstein IPG im¬ 
age (Figure 4 of Paper I) demonstrates that most of the 
features seen here were already present in the IPG image. 
The peculiar offset between the X-ray and the galaxy dis¬ 
tribution maxima noticed before - about 0.8 arcmin - is 
also evident. As discussed in Paper I, this offset is a con¬ 
sequence of the displacement of the X-ray isophotes in the 
direction of a concentration of galaxies to the northwest 
(a ~ 10*'16™48®, 6 ~ —10°38') which is centered on the 
brightest cluster galaxy. This galaxy group however has 
no associated X-ray excess (see also Fig|^, below). 

Given the moderate spatial resolution of the MEGS 
and LEGS detectors (FWHM 1.4 arcmin at 6 keV and 
3.6 arcmin at 2 keV, respectively), we have defined sub- 
regions in the field of view in order to perform separate 
spectral analysis. The X-ray data in these regions were 
extracted with XSELECT 2, from the HEAsoet package]^ 

We have defined 6 concentric annuli, each 1'28" 
wide (11 pixels), centered at the X-ray maximum (a = 
10^17™21®, 5 = —10°41'32"). These regions are employed 
to derive the radial prohles of various physical quantities. 
Notice that the last annulus extends over the detector re¬ 
gion affected by the MEGS “strongback”. 

We have also defined 4 regions corresponding to the 
NW, NE, SW, and SE quadrants (excluding the central 
region, see Sect. ||for details). Thi s subdivision in sectors 
i s very much like the one used by De Grand! fc Molendi 
( 1999|) . These regions, which are depicted in Fig. below, 
were chosen based on the galaxy distribution and their 
velocity distribution (c.f. Paper I). Finally, we have de¬ 
fined a circular region of 2 arcmin (15 pixels) radius cen¬ 
tered at the peak of the SE X-ray excess (a = 10^17™37®, 
5 = -10°46'02", see Fig. 0). 

The X-ray spectrum of each region was modelled as be¬ 
ing produced by a single temperature plasma and we em¬ 


ployed the MEKAL (Kaastra & Mewe 1993; Liedahl et al 


1995) model. For the whole cluster data, we also used the 
VMEKAL model which allows the individual abundances 
for the elements to be fitted independently. However, the 
abundance of some metals cannot be constrained by the 
spectral fitting, therefore we have only used the following 
elements: Si, S, Ar, Ga, Fe, and Ni. The others metals were 
fixed to O.SZq. 

The spectral fits were done usi ng XSPEC vll.O . 
Following the BeppoSAX “Gookbook” ( Fiore et al. 1999| ), 
we have selected the data in the energy range [0.12-4.0] 
keV for the LEGS, and [1.65-10.5] keV for the MEGS. 


® A service of the Laboratory for High Energy Astrophysics 
(LHEA) at NASA/ GSFC and the High Energy Astrophysics 
Division of the Smithsonian Astrophysical Observatory (SAO). 
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10hl8"i00s 17^45® 17^30® n^OO^ 16^45® 


Right Ascension (J2000) 


Fig. 1. X-ray isocontours obtained with BeppoSAX MEGS 2 and 3 superposed on a red DSS image (shown in grey 
scale). The larger dashed circle indicates the area used in our analysis; the smaller dashed circle indicates the excess X- 
ray emission tentatively associated with an active galactic nucleus (see text). The contours are spaced logarithmically. 
Notice the offset between the X-ray maximum and the centre of the galaxy distribution [a = 10^17™26.4®, 5 = 
—10°41'27", J2000, marked with an “X”, see also Fig. ^ below). 


The photoelectric absorption - mainly due to neutral hy- 
dr ogen - was computed using the cross-s ections given 
by Balucinska-Church fc McCammon ( 1992| ), available in 
XSPEC. Figure ^ shows four examples of spectral fits. 

For the last annulus, we have also done the spectral 
fits using only MFCS data in the interval [5.0-10.5] keV, 
in order to check the masking effect of the “strongback”. 
We verified that the difference between the best fit values 
obtained either in the full energy range [1.65-10.5] keV or 
in the [5.0-10.5] keV range is well within the error bars 
(cf. below, Sect. 3.2). 

Background spectra were constructed from the event 
files (produced by the BeppoSAX Science Data Centre) 
of the standard background data, which represent an as¬ 
semblage of several empty fields devoid of any notice¬ 
able source. For the MFCS, we have used the combined 
background for detectors 2 and 3 available from the 
BeppoSAX web site (under /pub/sax/cal/bgd/98_ll). 
For the LEGS, we have used the background event 
file described by Parmar et al. (1999) and available at 
/pub/sax/cal/bgd/99_12. The background spectra were 
extracted at the same regions (in detector coordinates) 
defined above and used in the spectral fits. 


3. Results 

3.1. Flux and luminosity 

Table summarizes the fluxes and luminosities obtained 
for the whole usable field and the central region. The bolo- 
metric X-ray luminosity implies Tx = 5.0 ±0.5 
on the Lx-Tx relation for clusters obtained by 

( W ). 


keV based 
Xue & Wu 


Table 1. Non-absorbed fluxes and luminosities in different 
energy bands. The energy band limits are given in keV. 


Region 

flux 

flux 

Lx 

Lx 

Lx 


[0.1-2.4] 

[0.5-4.5] 

[0.1-2.4] 

[0.5-4.5] 

bolom. 

All field 

27.8 

30.0 

42.2 

46.3 

83.1 

Centre 

5.13 

5.69 

7.78 

8.75 

16.3 


Note: flux units are 10 ^^erg cm 
Luminosity units are 10"*® erg s“^. 


As suggested in Paper I, the SE X-ray emission excess 
visible in Fig. may be due to an active nucleus excited by 
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channel energy [keV] 

Fig. 2. Four examples of the spectral fit. The left spectrum 
in each panel is from the LEGS, the right one from the 
combined MEGS 2 and 3. The lines passing through the 
data points correspond to the best fit absorbed MEKAL 
model. 


the merger of the pair of interacting galaxies. In order to 
verify this hypothesis, we have extracted a spectrum in a 
2 arcmin circle comprising this feature and tested different 
spectral models (see Table ||) . 


Table 2. Spectral fits of the south-east X-ray excess. Nh 
is the hydrogen column density in units of 10^°cm“^, kT 
is the temperature in keV, Z is the metallicity, a is the 
spectral index and “dof” is the degrees of freedom. Errors 
are at the Icr confidence level. 


Model 

Nh 

kT 

Z/Zq 

a 

xVdof 

POW 

< 4.1 

— 

— 

1.5 ±0.1 

80.7/89 

MEKAL 

< 2.5 

15.2 ±7.5 

0.35* 

— 

80.0/89 

MEKAL 

6.05* 

14.2 ±6.6 

0.35* 

— 

81.1/90 

MEKAL 

< 50 

4.5* 

0.55 ±0.29 

— 

94.1/89 

MEKAL 

6.05* 

4.5* 

0.58 ±0.33 

— 

95.3/90 


* Value fixed. 


We have tested a power-law (pow) and a thermal 
plasma (mekal) models. A composite model, power-law 
plus mekal (i.e., a non-thermal source superposed on the 
cluster thermal emission) could not be fitted unambigu¬ 
ously - fits with different initial parameters ended with dif¬ 
ferent minimum in parameter space. Therefore, we con¬ 
sidered two opposite scenarios: I) the X-ray excess comes 
from a non-thermal source and the emission inside a circle 
of 2 arcmin will be dominated by a power-law spectrum; 
II) the merging galaxies have no AGN and the emission is 
mostly of thermal origin. 

We found that the best models, both statistically 
equivalent, are the power-law and the thermal plasma en¬ 
dowed with the unreasonable high temperature of 14 keV 
or more. Fixing the plasma temperature at 4.5 keV (i.e., 
about the cluster mean temperature, see below), whether 
or not keeping fixed metallicities and Nh, produces poorer 
fits. We are thus led to suggest that the X-ray source is of 
non-thermal origin, probably due to the nuclear activity 
of one or both of the components of the interacting pair of 
galaxies optically associated to it. This conclusion should 
however be regarded with caution given the poor quality 
of the spectrum in this region, with net counts rates of 
5.1 X 10“^ count/s for the MEGS and 2.3 x 10“^count/s 
for the LEGS. 


3.2. Rddial profiles 


We have computed the radial profiles for several astro- 
physical quantities using the circular annuli defined in 
Section ^ with the combined MEGS and LEGS data. Since 
the SE X-ray excess may be of non-thermal origin, we have 
done the spectral fits excluding and including this region. 
The results are virtually the same, because the contribu¬ 
tion of the SE excess to the X-ray counts inside the annuli 
is not significant. The X-ray SE excess was therefore ex¬ 
cluded from the data prior to the spectral fits presented 
here. 


Figure]^ shows the resulting temperature profil^ 

The temperature profile was also determined using 
only MEGS da ta and fixing Nh = 5.42 x 10^°cm“^ (the 
galactic value, ( Dickey fc Lockman 1990 ). This was done 
for two reasons: first, to check the robustness of the fits, 
whether or not there was any spurious effect introduced 
by the use of two different detectors. 

The second reason is that when both the temperature 
and hydrogen column density are taken as free parame¬ 
ters, and depending on the spectral range being adjusted, 
the spectral fitting procedure may produce spurious cor- 
relations between best-fit values of kT and Nh [see, e.g., 
Pislar et al. ( 1997 )]. 


^ Actually, it is the emission weighted temperature that is ob¬ 
tained in a spectral fit. However, for physical realistic gas mod¬ 
els, the difference between the actual gas (electron) tempera¬ 
ture and the measured emission weighted temperature should 
be negligible and well w ithin the observational uncertainties 
(Komatsu & Seljak 2001). 
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ideal gas, P = nkT {n is the number density of particles), 
therefore the temperature and density are related by: 



0 

MECS+LECS 


0 

MECS full only 

- 

□ 

MECS 5-10+LECS 



-y= 1.16 


kT < 


, 7-1 


( 2 ) 



For the gas density profile, we use the well known /3- 
model (Cavaliere & Fusco-Femiano 1976): 


i(r) = no 


1 + 


-3/3/2 


( 3 ) 


4 6 

R (arcmin) 

Fig. 3. Temperature profile using either both LEGS and 
MEGS or only MEGS data (which is more appropriate 
for cluster temperatures T ^ 3 keV). For the last annulus 
there is an extra data point (square) where the spectral ht 
was done using only the MEGS [5.0-10.5] keV band. The 
error bars are at the Icr confidence level. The horizontal 
lines show the mean temperature (full line) and the la 
error bar (dashed) measured using all data within 8 arcmin 
[kT = 4A6't.Q'\‘^). The curve (full line) shows the best fit 
of a polytropic model (see text). 


where no, Tc and /3 are the central density, the core radius, 
and the shape parameter. We adopt the values determined 
in Paper I: /3 = 0.66 ± 0.44 and = (260 ± 20)/i^Q^kpc 
or 2.8 ± 0.2 arcmin. For our considerations we have used 
the Einstein IPG results because this instrument provides 
a higher spatial resolution (about 25 arcsec) compared to 
the BeppoSAX MEGS (about 1.5 arcmin). The tempera¬ 
ture profile is then written as: 


T(r) = To 


, r 

1 + ( - 

Tn 


n -3/3(7-1)/2 


( 4 ) 


where Tq is the central temperature. 

A least-square fit of Eq. (^) yields 7 = 1.15 ± 0.05 and 
To = 4.93 ± 0.23 with = 9.3 for 4 degrees of freedom. 


This value agrees well with those predicted by Komatsr 


& Seljak (2001) and found by Finoguenov et al 


( 2001 b). 


As it can be seen from Fig. the resulting temperature 
profiles, that is, for MEGS data only and /Vh fixed at the 
galactic value and for the combined MEGS and LEGS data 
and /Vh free, are almost indistinguishable. This demon¬ 
strates the internal consistency of the modelled temper¬ 
atures and a spurious correlation between fitted kT and 
/Vh values is not present. 

An additional fit for the outermost annulus was done 
using only the [5.0-10.5] keV MEGS data (but still using 
the full LEGS). This was done to verify whether or not 
the MEGS “strongback” has an important masking effect 
for low energy photons. It turns out that the effect on the 
temperature determination is within the error bars (which 
increases compared to the fit using the full MEGS range). 

A simple parametric form of the IGM temperature can 
be obtained using the polytropic equation of state. A poly¬ 
tropic model for the IGM is motivated by the results of 
hydrodynamical simulations (e.g., Suto et al. 1998|), the- 


& Dor^ 2002 |), and observed temperature gradients (e.g.. 

Markevitch et al 

1999 

Finoguenov et al. 

2001 a, , but see 

also Irwin et al.|| 

1999; 

rwin & Bregman 

2000 ; 

De Grandi 

& Molmdi'pOO^ 

for a different observational result - see 


discussion below). 

The polytropic gas is defined by: 

P (X p'^ , 


( 1 ) 


where P and p are the pressure and density of the gas, 
respectively, and 7 is the gas polytropic index. For an 


The fit, shown in Fig. is poor because of the very steep 
temperature gradient beyond 5 arcmin. In their study of 
a sample of 21 nearby clusters observed with BeppoSAX, 
De Grandi & Molendi ( 2002| ) hnd that the temperature 
profiles of clusters are characterized by isothermal cores 
extending to about one-hfth of the cluster virial radius, 
followed by steep temperature gradients which are steep¬ 
est for the non-cooling flow clusters. These authors also 
find poor fits when trying to adjust a polytropic law to 
the overall temperature profiles. This situation is much 
similar to that observed in Abell 970 which, when com¬ 
pared to their sample of non-c ooli ng flow clusters, presents 
a steeper gradient. In section 4.2 we use both the isother¬ 
mal and the polytropic model to estimate the total cluster 
mass. 

We may check whether the last three radial points still 
correspond to a stable state or whether the temperature 
gradient is too steep for stability. If we take into account 
only the three external points, i.e., 5 < R < 8 arcmin, 
the polytropic index increases to 7 = 1.7 ± 0.2, which is 
slightly above the ideal gas value, 5/3, suggesting that the 
outer region may be barely in stable equilibrium. 

The MEKAL plasma model spectral fit also gives the 
IGM metallicity. Figure ^ shows the metal abundance ra¬ 
dial profile. In order to check the ht robustness, we also 
did the prohle determination using only the MEGS data, 
besides the hts using both detectors. Both prohles are very 
similar, well within the error bars. 

As before, we made an additional ht for the last an¬ 
nulus using only photons with E > 5keV. Although the 
metallicity value obtained is close to the spectral ht with 
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Fig. 4. Metallidty profile using either both LEGS and 
MEGS or only MEGS data (the metallidty is determined 
mainly by the iron complex at ~ 6.8 keV). The square 
symbol has the same meaning as in the previous figure. 
The error bars are la confidence level. The horizontal lines 
show the mean metallicity measured using all data inside 
8 arcmin (Z = O.Slto.o^.^©). 


the full MEGS range, the error bars are so big that, in 
fact, the metallicity should be considered undetermined 
at this last annulus. 

The metallicity profile shown in Fig. ^ has no appar¬ 
ent gradient, being consistent with the mean value found 
for the whole cluster (which is a typical value for a clus¬ 
ter of such richness and temperature). Gontrary to some 
other clusters, e.g., Abell 496 (Dupke & White 2000, 


us¬ 


ing ASGA data) or Abell 85 (Lima Neto et al. 2001, using 
BeppoSAX data), this cluster does not show an increase of 
metal abu ndance in the cent r e (def ined by the X-ray peak 
emission). Irwin fc Bregman ( 2001 ), also using BeppoSAX 
data, find that clusters generally show a negative metal¬ 
licity gradient, but this gradient is usually smaller or even 
disappears in non-cooling flow clusters. 

Figure ^ shows the radial profile of hydrogen column 
density. For comparison we also show the galactic value 
at the position of Abell 970. Here we observe a gradient, 
with significantly higher hydrogen column density in the 
inner 2 arcmin. This excess is significant at a 2-a level, 
ignoring the uncertainty of the galactic iVu measure. 

This level of significance may indicate that we are de¬ 
tecting an X-ray absorption excess compared with the 
expected absorption predicted by the observed galactic 
HI column density. The origin of such an excess is still 
controversial (Allen 2000, and references therein), the 
most promising hypotheses being very cool molecular gas 
and/or dust grains. Moreover, the amount of this cool 
material depends on the details of the spectral modelling 
( McGarthy West fc Welch 2002| ). 

Figure || shows that iVn decreases outwards, reaching 
the galactic value. However, a caveat must be introduced 
here: (I) there may be variations of the galactic A^h that 
are unresolved by the galactic survey at the cluster scale; 


Fig. 5. Hydrogen column density profile using LEGS and 
MEGS data. The error bars are Icr confidence level. Filled 
diamonds correspond to fits with LEGS counts restricted 
to the range [1-4] keV (but all MEGS). The open square 
corresponds to the last annulus fit in which MEGS counts 
were restricted to the range [5-10.5] keV (see text). The 
horizontal full line shows the mean hydrogen column den¬ 
sity measured using all data within 8 arcmin; the hori¬ 
zontal dashed lines correspond to Icr error bars (A^h = 
6 . 051/097 ^ 10^°cm“^). The horizontal dotted line indi¬ 
cates the mean galactic value, A^h = 5.42 x 10^°cm“^. 


( 11 ) the A^h determination relies almost exclusively on the 
LEGS data (photon energy ^1.5 keV). The error bars in 
this case may be underestimated given the low count rate, 
specially at the last bin. 

Moreover, there are a number of issues associated with 
lowest energy photons detected by the LEGS: the spatial 
resolution is poor (« 5.14/'\/£'kev arcmin) and there is a 
disagreement between the measured (on-gr ound) and the¬ 


oretical effective area for E ^ 0.3 keV (see Parmar et al 


19971) . Selecting only LEGS counts with E > 1.0 keV (the 
spatial resolution is ~ 5 arcmin at if = 1 keV) can mini¬ 
mize these effects. Figure ^ shows the additional A^h de¬ 
termination with such a selection, for the centre and the 
first two rings. The error bars are significantly larger, but 
the central excess in modelled neutral hydrogen column 
density is still present at 2 ct significance. 


4. Cluster mass 

4.1. Gas mass 

In order to determine the cluster gas mass we assume that 
the gas distribution can be described by a spherically sym¬ 
metric 6 efa-model profile with values of core radius, r^, 
and shape parameter, (3 as determined in Paper 1. With 
the good quality spectrum we now dispose, it is possible 
to estimate the central electronic density, uq, as follows. 

We integrate the bremsstrahlung emissivity along the 
line-of-sight, in the energy bands given in Table ^ and 
obtain uq. Taking into account the error bars in and 
/3, we obtain no = (3.7 ± 0.6)10“^cm“^. Notice that the 
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shape of the temperature profile has only a minor effect 
(less than 5%) in the determination of ng. 

The integrated gas mass profile (or growth curve) has 
been obtained by integration of Eq. (^ and is shown below 
in Fig. 


4.2. Dynamical mass 


With the density and temperature profiles, we can esti¬ 
mate the total cluster mass. Supposing the ICM is in hy¬ 
drostatic equilibrium, the total “X-ray” dynamical mass 
(i.e., the dynamical mass determined from an X-ray ob¬ 
servation, not to be confused with the mass of the X-ray 
emitting gas) inside r is given by: 


M{r) 


kT 

G/iTOp 


/d In p{r) 
y dlnr 


dlnr(r) \ 
dlnr / 


(5) 


where fi rup is the mean mass per particle (/i = 0.59 for a 
fully ionized primordial gas) and nip is the proton mass. 
We assume that matter is distributed spherically. 

Adopting the /3-model to describe the intracluster gas 
density and the polytropic equation of state, the dynami¬ 
cal mass results in: 


M{r) 


SkToPyrc 

Gprup 





( 6 ) 


When 7 = 1 we have the usual expression for the 
isothermal sphere. Figure shows the mass profile com¬ 
puted both with an isothermal temperature and poly¬ 
tropic model. 

The dynamical masses estimated here are about 15% 
to 25% higher than those given in Paper I where a 3.3 keV 
isothermal model was assumed. 

We can now compute the radius corresponding to the 
density ratio p{rs)/pc{z) (see Appendix). For S = 200, 
we have the usual r 2 oo that is ofte n used as an approx - 
imation to the virial radius (e.g., Lacey fc Cole 1993|) . 
Assuming either an isothermal or polytropic profile, we 
obtain r 2 oo = 2.3/i^g^Mpc and r 2 gg = 1.9/i^g^Mpc, re¬ 
spectively. Thus, i?image) the maximum radius analysed, 
is 0.33-0.4 r 2 gg, depending on the adopted temperature 
profile. 

From the dynamical and gas mass profiles, we can de¬ 
rive the baryon fraction profile, /(r) = MbaryonZ-ZLltotai, 
assuming that the bulk of the baryons are in the cluster 
ICM (i.e., Mbaryon ~ Mgas). The neglect of the galactic 
contribution (about ^ 5 to 10% of the total mass, see 
Paper I), underestimates / by just a few percent. 

The bottom panel of Fig. gives the estimated baryon 
fraction profile. For reference, we show r' 25 gg and rgggg, 
the radii within the mean cluster density is 2500 and 5000 
times the value of the critical density of the Universe. 

The gas mass fraction obtained here agrees with the 
typical values obtained by White fc Fabian] (1995) (at R = 
1 Mpc) and by Allen et ^ (2002) (at R = ?’ 25 gg). However, 
the shape of the /(r) is much flatter than the profiles 
computed by Allen et ah. A possible explanation might 



0 2 4 6 8 10 

R (arcmin) 


Fig. 6. Top: X-ray dynamical mass prohle using an 
isothermal model, with the mean cluster temperature 
(solid line), and polytropic model, with central temper¬ 
ature Tg = 5.0 keV and 7 = 1.15 (dashed curve). The 
gas mass is also shown (dotted curve). Bottom: Baryon 
fraction given by the dynamical and gas mass ratio. It is 
shown r 25 gg and r^ooo: where the density contrast, p{r)/pc, 
is 2500 and 5000, respectively (see Appendix). i?image is 
the maximum radius analysed. 


be the adopted model for t he mass: wh ile we assume a 
/3-mode l for the gas profile , Allen et al. adopt the NFW 
profile ( Navarro et al 1997 ) for the total mass. 

Finally we must stress that although the new estimates 
of the X-ray dynamical mass presented here are higher 
than those given in Paper I, they are still a factor, ^ 5- 
10, below the virial mass estimates (VME) obtained by the 
application of the Virial Theorem to the galaxies. As dis¬ 
cussed in Paper I, in the case of equilibrium, the dynamical 
mass by the VME should not be overestimated by more 
than 20%. We are thus lead to interpret the large discrep¬ 
ancies between the virial theorem and the X-ray dynami¬ 
cal mass estimates as evidence of a non-equilibrium state 
of cluster. Indeed, in Paper I we had already remarked the 
incidence of substructures in the galaxy projected distri¬ 
bution and of a large scale gradient of their line-of-sight 
velocities, which are strong signs that the cluster is out of 
equilibrium. We will discuss more in the following Section. 


5. Cluster merging 

Figure ^ shows the temperature and metallicity measured 
in four quadrants: NW, NE, SW, and SE region. 

As it can be seen from this figure, there is a hint of a 
temperature gradient in the North-South direction. Also, 
there is a possible metallicity gradient in the NE-SW di- 
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Right Ascension (J2000) 


Fig. 7. Temperature and metallicity in the regions cen¬ 
tral, NW, NE, SW, and SE, superposed on the combined 
MEGS 2 and 3 adaptively smoothed image. The solid lines 
mark the cluster regions for which we carried out the spec¬ 
troscopic analysis. The central 1'28" region corresponds to 
the innermost point in Figs. The isocontours (spaced 
logarithmically) of the projected density map of galaxies 
kinematically bound to the cluster and with B < 19, are 
superposed to the X-ray image (cf. Fig 1 in Paper I). 


rection. However, within their error bars, those values are 
also compatible with an isotropic distribution. In Paper I, 
it has been shown that there is a large-scale gradient of 
the local mean velocities of the cluster galaxies. The gen¬ 
eral direction of this velocity gradient is approximately 
the same as the temperature gradient hinted in Figure |^. 

Moreover, we observe that the irregular galaxy distri¬ 
bution does not coincide with the X-ray emission map. 
There is also a discrepancy between the dynamical mass 
determined by the virial theorem and X-ray gas emission 
(see Sect. These features are recognized as signs of past 
mergers between sub-clusters. We may infer that, in this 
case, the merger is not so recent, since the X-ray distribu¬ 
tion already shows a fairly axial symmetric distribution. 

From another perspective, as discussed in Paper I and 
shown in Fig. the galaxy distribution presents a sub¬ 
structure located at a projected distance of 760/i^o^kpc 
(~ 8 arcmin) NW of the centre, with a velocity dispersion 
of 380t7o°km s“^ inside a circular radius of 285 /i^Q^kpc 
(^ 3 arcmin). This was interpreted as a group falling into 
the core of Abell 970. Taking into account the mean veloc¬ 
ities of the central region, 17600± 118 km/s, and NW sub¬ 
structure, 17834± 135, (cf. Paper I), then we can roughly 
assume that the substructure may be falling with a veloc¬ 
ity in the range 50-350 km/s. Consequently, with such a 
crude estimate, we conclude that a merger will occur in 
the next 2-6 Gyr. 


Another possibility is that this substructure has al¬ 
ready passed near the cluster centre, disturbing the over¬ 
all cluster equilibrium, and will come back again to the 
cluster central regions, before being destroyed. 

Therefore, the galaxy distribution (spatial and veloc¬ 
ity) together with the X-ray isophote off-centre and pos¬ 
sible temperature gradient suggest that Abell 970 is be¬ 
tween two merger events: the first one that is responsible 
for the offset between the X-ray and galaxy distributions; 
the second one will occur when the NW substructure falls 
into the central part of the cluster. 


6. Metal abundances 


The pattern observed in the heavy-element abundances 
of the intracluster medium provides important clues on 
the types of galaxies and stars producing these elements 
and the mechanisms responsible for their removal from the 
interstellar medium of their parent galaxies. 


Using the whole 8 arcmin field (0.33-0.4 r 2 oo, see 
Sect. 4^), we fitted a plasma model with the following 
metal abundances as free parameters: Si, S, Ar, Ca, Fe, 
and Ni. The following element abundances were fixed to 
0.3^0: He, C, N, O, Ne, Na, Mg, and Al. 

The result is summarized in Figure which shows 
the over-abundance of a-elements relative to iron. Both 
the silicon and sulphur abundances are consistent with 


t he typ ical values found recently by Finoguenov et al 
( 2001 b|) and are slight ly above the value determined by 


Fukazawa et al. ( |1998D for a cluster of ~ 4.5 keV (which 
is ~ 0.6-0. 8 ^ 0 , cf. Fig. 2 of their paper). These authors 
have shown that the Si abundance correlates with the gas 
temperature, probably due to a decreasing contribution 
from Type H supernovae (SNe) in poorer clusters, where 
the metals produced by this type of SN could escape from 
the potential well. 


The over-density of a-elements displayed in Fig. || is 
characteristic of the Type H SNe yield. This suggests that 
the intracluster gas enrichment was caused by ejected gas 
that was synthesized in massive stars, the progenitors of 
Type H SNe. Present day elliptical galaxies show a very 
low rate of Type H SNe, implying that the enrichment of 
the intracluster gas probably occurred early in the history 
of the cluster. 


However, we cannot discard other possible enrichment 
mechanisms. In fact, as shown in Figure |^, the abundance 
ratio Ni/Fe is higher than expected for Type H SNe, but 
it is in agreement with the expected Type la SNe yield. 
This finding supports models where part of the metals in 
the ICM were produced by Type la SNe and which were 
either blown from the galaxies by SNe explosions or swept 
by ram-pressure (e.g., Dupke fc Whit^ ( p000| )). However, 
the Ni abundance is uncertain. Moreover, the abundance 
ratio Si/Fe (with its smaller error bars) is incompatible 
with enrichment by Type la SNe. 
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Fig. 8. Top: Individual metallicities of Fe, Ni and a- 
elements obtained with an absorbed VMEKAL model. The 
solid horizontal line is the mean metallicity obtained with 
an absorbed MEKAL model; the dashed lines correspond 
to la confidence level, except for argon which is an upper 
limit. Bottom: abundance ratios between metal and iron. 
Also shown are the predicted approximate yields of Type 
II SNe (diamonds) and Type la SNe (crosses), based on 
Dupke fc Whit^ ( pOOO ). 


the nearby NW substructure will fall (again?) into the 
core of Abell 970. 

— We have derived an abundance ratio between a- 
elements/iron greater than 1 and the iron metallicity is 
lower than the mean metallicity. The over abundance 
of a-elements may be an indication of an early enrich¬ 
ment of the ICM by Type II supernovae in elliptical 
galaxies. Part or the metals in the ICM could be pro¬ 
duced by Type la SNe, but the evidence from the Ni/Fe 
ratio is not strong. 

— We have found a central X-ray absorption excess com¬ 
pared to the absorption expected from the mean galac¬ 
tic value of iVn, although the statistical error bars are 
large. The central iVn excess could be associated to the 
central galaxies and/or to cooled intra-cluster material 
(either in molecular form and/or dust grains). The lat¬ 
ter may be preferable since the core of the cluster is 
dominated by early-type galaxies. However, this result 
probably depends on the spectr al modelling (see, e.g., 
[McCarthy, West fc Weldj ^002 ). 
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7. Conclusions 

We have presented here a new X-ray observation of Abell 
970. Below we summarize our main results. 


Appendix A: Derivation of rs 

In Paper I, we give a formula for rs, which is defined by: 


— Although Abell 970 presents rather symmetric X-ray 
isophotes (except for an excess emission that is prob¬ 
ably due to an active galactic nucleus, triggered by a 
galaxy merger), the peak of the X-ray distribution does 
not coincide with the centre of the galaxy distribution. 

— The temperature profile presents a strong radial gra¬ 
dient at ^ 6 arcmin, which is barely compatible with 
adiabatic stable equilibrium. The metal abundances 
are approximately constant up to 8 arcmin. There are 
indications of non radial temperature and metallicity 
gradients (but with conhdence levels less than Icr), 
notably in the NW-SE direction, the same direction 
where a gradient in velocity is observed (cf. Paper I). 

— With the two findings above, we arrive at two possible 
scenarios for the recent history of Abell 970: either it 
has suffered a recent merger with a subcluster, or the 
NW substructure (detected only from galaxy counts) 
has recently passed through the centre of Abell 970. 
In both cases, the dynamical stability was disrupted, 
but the cluster is now approaching virial equilibrium 
again. The observed galaxy distribution suggests that 
this equilibrium will be broken in a few Gigayears when 


c _ P{rs) 

0 = —TT , 

Pc{z) 


(A.l) 


where p{rs) = Mdyn(r5)/(47rr|/3) is the mean density of 
the cluster inside rs and pdz) is the critical density of the 
Universe at the redshift of the cluster. 

We take the dynamical mass corresponding to a /3- 
model ICM and polytropic temperature prohle, Eq. (^, 
and the critical density given by: 


Pc(^) = ^^^; H{z)=Hof{z,nM,n^), (A.2) 

with 

P{z,nM,^A)= /.O', 

Ha + Hm(1 + z)^ - (Ha + Hm - 1)(1 + z)^ . ^ ’ 


Combining these equations and solving for rs we obtain: 



2.3x10® /3jTo 
ShlQp{z,nM,^A) prl 


1 /[ 1 + 3 / 3 ( 7 - 1 )/ 2 ] 


(A.4) 
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where the cluster central temperature, Tq, is given in keV 
and the c ore radius, Vc, in kpc. For <5 ^ 500, rs ^ fc, and 
Eq. (A.4) is simplified: 


rs_ 

Ta 


2.3x10® 


5hl^P{z,Q.M,^K) 


(A.5) 


When 7 = 1, the isothermal case is obtained and Tq is the 
mean temperature. 

For z 0.1 we can use the approximation: 
f^{z, Umi Ha) — 1 + (2 — 2r2A + ^m) z . (^^-6) 


References 

Allen S.W., 2000, MNRAS 315, 269 
Allen S.W., Schmidt R.W., Fabian A.C., 2002, MNRAS 
334, Lll 

Balucinska-Church M., McCammon D., 1992, ApJ 400, 
699 

Boella G., Butler R.C., Perola G.C., Piro L., Scarsi L., 
Bleeker J.A.M., 1997a, A&AS 122, 299 
Boella G., Ghiappetti L., Gonti G., Gusumano G., del 
Sordo S., La Rosa G., Maccarone M.G., Mineo T., 
Molendi S., Re S., Sacco B., Tripiciano M., 1997b, 
A&AS 122, 327 

Gavaliere A., Fusco-Femiano R, 1976, A&A 49, 137 
Gavaliere A., Mend N., Tozzi P., 1999, MNRAS 308, 599 
Gusumano G., Mineo T., 1998, “On-flight calibration of 
the MEGS vignetting” 

De Grandi S., Molendi S., 1999, ApJ 527, L25 
De Grandi S., Molendi S., 2002, ApJ 567, 163 
Dickey, Lockman 1990, Ann. Rev. Ast. Astr. 28, 215 
Dos Santos S., Dore O., 2002, A&A 383, 450 
Dupke R.A., White R.E., III, 2000, ApJ 537, 123 
Einasto M., Tago E., Jaaniste J., Einasto J., Andernach 
H., 1997, A&AS 123, 119 

Evrard A.E., MacFarland T.J., Couchman H.M.P., 
Colberg J.M., Yoshida N., White S.D.M., Jenkins A.R., 
Erenk C.S., Pearce F.R., Peacock J.A., Thomas P.A., 
2002, ApJ 573 

Finoguenov A., Reiprich T.H., Bohringer H., 2001a, A&A 
368, 749 

Finoguenov A., Arnaud M., David L.P., 2001b, ApJ 555, 
191 

Fiore F., Guaianazzi M., Grandi P., 1999, “Gookbook 
for BeppoSAX NFI spectral analysis”, version 1.2, 
WWW.sdc.asi.it/pub/sax/doc/software_docs/saxabc_vl. 2 
ps.gz 

Fukazawa Y., Makishima K., Tamura T., Ezawa H., 
Haiguang X., Ikebe Y., Kikuchi K., Ohashi T., 1998, 
PASJ 50, 187 

Irwin J.A., Bregman J.N., Evrard A.E., 1999, ApJ 519, 
518 

Irwin J.A., Bregman J.N., 2000, ApJ 538, 543 
Irwin J.A., Bregman J.N., 2001, ApJ 546, 150 
Jones G., Eorman W., 1999, ApJ 51, 65 
Kaastra J.S., Mewe R., 1993, A&AS 97, 443 


Kalinkov M., Valtchanov L, Kuneva I., 1998, A&A 331, 
838 

Komatsu E., Seljak U., 2001, MNRAS 327, 1353 
Lacey G., Gole S., 1993, MNRAS 262, 627 
Lea S.M., Silk J., Kellogg E., Murray S., 1973, ApJ 184, 
105L 

Liedahl D.A., Osterheld A.L., Goldstein W.H., 1995, ApJ 
438, L115 

Lima Neto G.B., Pislar V., Bagchi J., 2001, A&A 368, 440 
Loken G., Melott A.L., Miller C.J., 1999, ApJ 520, L5 
Markevitch M., Vikhlinin A., Forman W.R., Sarazn G.L., 
1999, ApJ 527, 545 

McGarthy LG., West M.J., Welch G.A., 2002, ApJ 567, 
762 

Navarro J.F., Frenk G.S., White S.D.M., 1997, ApJ 490, 
493 

Parmar A.N., Martin D.D.E., Bavdaz M., Favata F., 
Kuulkers E., Vacanti G., hammers U., Peacock A., 
Taylor B.G., 1997, A&A Suppl. 122, 309 
Parmar A.N., Oosterbroek T., Orr A., et ah, 1999, A&A 
Suppl. 136, 407 

Pislar V., Durret F., Gerbal D., Lima Neto G.B., Slezak 
E., 1997, A&A 322, 53 

Schwope A., Basinger G., Lehmann L, Schwarz R., 
Brunner H., Neizvestny S., Ugryumov A., Balega Yu., 
Trumper J., Voges, W., 2000, AN 321, 1 
Sodre L., Jr., Proust D., Gapelato H.V., Lima Neto G.B., 
Guevas H., Quintana H., Fouque P., 2001, A&A 377, 
428 (Paper I) 

Suto Y., Sasaki S., Makino N., 1998, ApJ 509, 544 
White D.A., Jones G., Forman W., 1997, MNRAS 292, 
419 

White S.D.M., Rees M.J., 1978, MNRAS 183, 341 
White D.A., Fabian A.G., 1995, MNRAS 273, 72 
Xue Y.J., Wu X.P., 2000, ApJ 538, 65 






